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Evidence for crust cooling in the transiently accreting 
11-Hz X-ray pulsar in the globular cluster Terzan 5 
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ABSTRACT 

The temporal heating and subsequent cooling of the crusts of transiently accreting 
neutron stars carries unique information about their structure and a variety of nu- 
clear reaction processes. We report on a new Chandra Director's Discretionary Time 
observation of the globular cluster Terzan 5, aimed to monitor the transiently accret- 
ing 11-Hz X-ray pulsar IGR J17480-2446 after the cessation of its recent 10- week long 
accretion outburst. During the observation, which was performed ~ 125 days into 
quiescence, the source displays a thermal spectrum that fits to a neutron star atmo- 
sphere model with a temperature for an observer at infinity of kT°° ~ 92 eV. This is 
~ 10% lower than found ~ 75 days earlier, yet ~ 20% higher than the quiescent base 
level measured prior to the recent outburst. This can be interpreted as cooling of the 
accretion- heated neutron star crust, and implies that crust cooling is observable after 
short accretion episodes. Comparison with neutron star thermal evolution simulations 
indicates that substantial heat must be released at shallow depth inside the neutron 
star, which is not accounted for in current nuclear heating models. 

Key words: globular clusters: individual (Terzan 5) - X-rays: binaries - stars: neutron 
- pulsars: individual (CXOGClb J174804.8-244648, IGR J17480-2446) 



1 INTRODUCTION 

Neutron stars are the densest directly observable stellar 
objects in our Universe and constitute ideal astrophysical 
laboratories to study matter under extreme physical con- 
ditions. The outer layer of a neutron star, its crust, covers 
about one tenth of the total stellar radius and consists of 
ions, electrons and neutrons. The structure and composi- 
tion of the crust play an important role in the emission of 
gravitational w aves and the ev olution of th e neutron star's 
magnetic field (Brown & Bildstcn 1998; Ushomirsk v et al.l 
|2000J ; iHorowitz fc Kadaull2009l '). Furthermore, studying the 
crusts of neutron stars allows th e investigation of a vari- 
ety of nuclear reacti on processes i|Haensel fc Zdunikl I2008J ; 
IHorowitz et~al"1 120081 ). 

When residing in a low-mass X-ray binary (LMXB), a 
neutron star strips off the outer gaseous layers of a stellar 
companion through Roche-lobe overflow, and this matter is 
accreted onto the neutron star. Transient systems typically 
exhibit outbursts that last several weeks, after which accre- 
tion onto the neutron star ceases and the binary returns 
to a quiescent state. Usually it remains as such for many 
years before entering a new active phase. The accretion of 
matter compresses the neutron star crust, which induces 



e-mail: degenaar@uva.nl 



a chain of exothermic electron captures, neutron emissions 
and density-driven nuclear fusion reactions that take place 
at several h undreds meters depth (|Haensel fc Zduniklll990l . 
2008; Gu pta et al.l 12007). The energy released in these re- 
actions locally heats the crust, but is eventually conducted 
over the entire stellar body and radiated away via neutri- 
nos emitted from the dense core, and vi a thermal X-ray ra- 
diation from the ne utron star surface IjBrown et al.l Il998l ; 
iRutledge et al.l 2002). During quiescent episodes, the ther- 
mal X-ray radiation emerging from the hot neutron star can 
be detected with sensitive X-ray satellites. 

Prior to the onset of an accretion phase, a neutron star 
is nearly isothermal and the surface radiation the refore offers 
a pro be of the temperature of the stellar core (IColpi et al.l 
l200ll ). Shortly after an accretion outburst, however, the 
thermal X-ray emission tracks the temperat ure of the heate d 
crust (jUshomirskv fc Rutledgej|200ll ; IRutledge et al.l [2002). 
As the crust cools and restores thermal equilibrium with 
the core, the thermal X-ray emission decreases and set- 
tles at a quiescent b ase level set by the core temperature 
IjRutledge et al.ll2002r ). This thermal relaxation depends on 
the composition of the crust, which determines the conduc- 
tion of heat towards the surface and the core, and on the 
details of the nuclear reaction s, which set the magnitu de 
and depth of the heat sources (|Brown fc Cummindl2009l ). 

In the past decade, crust cooling has successfully been 
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monitored for four transiently ac c reting neutron stars : 
KS 1731-260 dWiinands etaflfeoOll; [Cackett et al.ll2010a[) , 
MXB 1659-29 IWiinands et all 12003: ICackett et alJl200a i. 
XT E J1701-462 llFridriksson et alj|201ll) and EX O 0748- 
676 JDegenaar et al.ll201ll ; iDiaz Trigo et al.ll201ll ). All four 
belong to a rare sub-group of transient LMXBs, the so-called 
quasi-persistent sources, which exhibit unusually long ac- 
cretion episodes extending up to decades. These served as 
prime targets for crust cooling studies, since the prolonged 
accretion ensures a substantially heated crust. Comparing 
the X-ray observations with thermal evolution simulations 
has provide d unique insight into the properties of neutro n 
star crusts (jShternin et al.ll2007l ; brown fc Cumming||2009h . 

An important conclusion drawn from these studies is 
that the crust efficiently conducts heat suc h that the ion 
lattice must have a h ighly ordered structure (jShternin et al.l 
120071 : brown fc Cum ming 2009). This poses problems for 
explaining the ignition of thermonuclear superbursts ob- 
served from some accreting neutron stars: these require a 
critical crust temperature that may be difficult to achieve 
when heat is efficiently c onducted (|Cooper fe Naravanl 
120051 ; ICumming et al.l [2006). As a possible solution, it 
has been hypothesized that accretion deposits more heat 
inside the n eutron star than is acc o unted for in cur- 
rent models J Brown fc Cummind 120091 : iGupta et all 120071 : 
borowitz et al.l 120081 ). New observations of neutron star 
crust cooling can confirm the existence, and constrain the 
magnit ude and location, of such p ossible additional sources 
of heat (JBrown fc Cummindl2009l ) . Unfortunately, X-ray bi- 
naries with long accretion phases are very rare: only a hand- 
ful of sources se rve as potential targets for future studies 
l|Wiinanddl200l ). 

IGR J17480-2446 (CXOGClb J174804.8-244648; J1748 
hereafter) is a transient LMXB located in the dense core of 
the globular cluster Terzan 5 (Fig. HI . It was initially discov- 
ered during Chandra observations of the cluster performed 
in 2003 and classified as a low- luminosity X-ray sou rce, pos- 
sibly a dormant X-ray binary iJHeinke et al.l l2006b). In 2010 
October, the s ource intensit y suddenly increased by orders 
of magnitude (JBordas et allboiol : boolev et ai]|20ich. The 
detec t ion of thermonucle ar X-ray bursts (jChenevez et al.l 
2010; Linares et alJ 20111) and coherent X -ray pulsations 
i Strohmaver et al.l2010l : IPapitto et alJ201ll ) unambiguously 
identified J 1748 as a neutron star LMXB. Timing analy- 
sis of the 11-Hz X-ray pulsations, representing the neutron 
star's spin period, revealed an orbital period of ~ 21 h and 
a ~ 0.4-1.5 M Q donor star (|Papitto et al.ll201lf) . 

In 2010 December, ~ 10 weeks after its onset, the accre- 
tion activity ceased and the b inary returned to quiescence 
l|Degenaar fc Wiinandsll2011al ). A Chandra observation ob- 
tained about 50 days later, in 2011 February, demonstrated 
that the quiescent thermal X-ray emission of J1748 was ele- 
vated above the base level measured from archival observa- 
tions taken in 2003 and 2009, by approximately a factor of 4 
l|Degenaar fc WiinandsH2011al lbr) . The associated rise in neu- 
tron star temperature (a factor of 1.4), can be interpreted as 
an accretion-heated crust, although two alternative explana- 
tions could be invoked (Degenaar & Wiinands 2011a). First, 
the neutron star possibly continues to accrete at a very slow 
rate, causing a variable quiescent flux. Second, the amount of 
hydrogen and helium left on the neutron star after an accre- 
tion phase determine the heat flux that is flowing from the 
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Figure 1. Chandra/ ACIS image of Terzan 5 obtained on 2011 
April 29-30. J1748 is located in the dense cluster core. 



stellar interior towards the surface. For a constant interior 
temperature, the quiescent thermal emission might therefore 
differ by a factor of a few due to changes in the envelope com - 
position after an intervening outburst (|Brown et al.l 120021 ). 
Both of these scenarios could also account for the enhanced 
surface temperature observed in 2011. 



2 OBSERVATIONS, ANALYSIS AND RESULTS 

2.1 A new Chandra observation of Terzan 5 

To further investigate the cause of the elevated thermal X- 
ray emission detected from J1748 after its 10- week outburst, 
we obtained a new Chandra Director's Discretionary Time 
observation of Terzan 5 on 2011 April 29-30 (obs ID 13252; 
Fig. [T]). This is ~ 125 days into quiescence and ~ 75 days 
since the previous Chandra observation. The cluster was po- 
sitioned on the back-illuminated S3 CCD at the nominal tar- 
get position and the observation was read out in the FAINT 
timed data mode. We used the CiAO software tools (v. 4.3) 
for the reduction and extraction of data products, following 
standard analysis procedures. There were no occurrences of 
background flares during the observation, so that the final 
net exposure time on our target was 39.5 ks. 

Using the task DMEXTRACT, we extracted source 
count rates and light curves from a circular region with 
a radius of 1", c entered at the b est-known position of 
IGR J 17480-2446 l|Poolev et alj|2010h . Corresponding back- 
ground events were collected from a circular region with a 
radius of 40" , positioned on a source-free part of the CCD 
located 1.4' west of the cluster core. J1748 is detected at a 
net count rate of (3.7 ± 0.3) x 10 -3 counts s _1 and a total 
of 146 source photons were collected during the observa- 
tion. Source and background X-ray spectra were obtained 
using the meta-task SPECEXTRACT, which also generates 
the appropriate ancillary resposonse files (arf) and redistri- 
bution matrix files (rmf). Using the FTOOL GRPPHA, the 
spectral data was grouped into bins with a minimum number 
of 20 photons. The resulting background-corrected spectrum 
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Figure 2. Chandra/ ACIS spectra of IGR J17480-2446 at four 
epochs. The 2011 data was obtained within a few months after 
the end of the 2010 October-December outburst and the 2009 
data 1.2 yr prior to the accretion activity. The solid lines indicate 
best-fits to the neutron star atmosphere model NSATMOS. 



was modeled in the energy range of 0.5-10.0 keV using the 
software package XSPEC (v. 12.6; lArnaudlll996r ). 



2.2 Quiescent spectra of IGR J 17480-2446 



Simila r to previous observations IjDegenaar fc Wiinandsl 
l2011al lbl). the X-ray spectrum of J1748 can be adequately fit 
by an absorbed thermal emission model and does not require 
the addition of any hard, non-thermal component (Fig. [2]). 
Within XSP EC, we use the neutr on star atmosphere model 
NSATMOS |Heinke et all 1200681 ) to fit the thermal emis- 
sion from the neutron star surface. In this model, we fix 
the mass and radius of the neutron star to M — 1.4 M© 
and R — 10 km, respectively, whereas the source distance is 
set to the current be st estimate for Terzan 5: D = 5.5 kpc 
lOrtolani et al.l 120071 ). The normalization was frozen at the 
recommended value of 1, which implies that the radiation is 
emerging from the entire stellar surface. Our only free pa- 
rameter for the NSATMOS model is thus the neutron star 
temperature. To take into account the hydrogen absorption 
along the line of sight, Nh, we include the PHABS model. 

To chart the evolution in thermal emission and neutron 
star temperature, we fit the spectrum of this new obser- 
vation together with Chandra data obtained on 2003 July 
13-14, 2009 July 15-16 and 2011 February 17. These had 
the same setup as the new observation and were treated fol- 
lowing the same reduction and analysis s teps used here (for 
details, see lDegenaar fe Wiinandsll2011al lbl). Following pre- 
vious work , we treat the 2003 and 2009 data as one single 
spectrum (JDegenaar fc Wiinandsl 1201 lbl ). We tied the hy- 
drogen column density between the data sets, i.e. this fit 
parameter is required to be the same amongst the obser- 
vations, whereas the neutron star temperature could vary 
freely for the different spectra. This approach resulted in a 
joint fit value of N H - (2.0 ±0.3) x 10 22 cm" 2 and a reduced 
chi-square of xt = 0.77 for 21 degrees of freedom. Fig. [5]dis- 
plays the X-ray spectra observed at the four epochs. 

For each data set we deduce the neutron star tempera- 



ture as seen by an observer at infinity, kT°° , and calculate 
the total unabsorbed model flux in the 0.5-10 keV energy 
range, F una b s . The model fits are then extrapolated to the 
0.01-100 keV energy range to obtain an estimate of the to- 
tal emitted bolometric flux, which is translated into a lu- 
minosity, Lboi, by assumi ng that the source is located at a 
distance of D — 5.5 kpc (|Ortolani et al.l 120071 ). The results 
of our spectral analysis are summarized in Table [TJ where 
all quoted errors refer to 90% confidence levels. The evolu- 
tion of the neutron star temperature as inferred from the 
spectral fits is plotted in Fig. [3] 



3 DISCUSSION 

We analysed a new Chandra observation of Terzan 5, ob- 
tained ~ 125 days after the end of the accretion outburst 
of J1748. The quiescent X-ray spectrum of the source is 
soft and can be adequately fit by a neutron star atmo- 
sphere model NSATMOS with a temperature of kT°° = 
92.1 ± 3.1 eV. This is nearly 10% lower than measured 75 
days earlier in 2011 February (kT°° = 100.9 ± 3.3 eV), yet 
over 20% higher than the quiescent base level of the source 
detected in 2003 and 2009 (fcT°° = 72.0 ±3.3 eV). The ther- 
mal X-ray emission and inferred neutron star temperature 
have thus decreased after the initial enhancement (Fig. [3} . 

To explain the elevated thermal emission observed in 
2011 in terms of a different envelope composition left after 
the 2010 accretion phase, the neutron star temperature in- 
ferred from our new observation should have been similar 
to the previous measurement, since the comp osition of the 
layer would not change in absence of accretion (|Brown et al.l 
|2002j). This is not consistent with our results. Continued ac- 
cretion is thought to reveal itself through stochastic vari- 
ability in the X-ray emission on both long (months/years) 
and short (seconds/days) time scales, and/or via the pres- 
ence of a hard, non-thermal emission component in the X- 
ray spectrum JRutledge et all l200ll ; ICackett et ail l2010bl ; 
iFridriksson et al.ll201ll ) . We find none of such features in the 
Chandra data of J1748, so there are no obvious indications 
that the neutron star is accreting at a low level[j 

On the other hand, the observed decrease in thermal 
X-ray emission and neutron star temperature is very remi- 
niscent of the crust cooling observed from the four X-ray bi- 
naries with long accretion phases. The data points of J1748 
(Fig. [3} are consistent with an exponential decay with an 
e-folding time of ~ 100 — 200 days or a power law with a 
decay index of order ~ 0.1. This is similar to the results ob- 
tained for the four quasi-persistent LM XBs (Cackct t et al.l 
120081 . l2010al ; iDegenaar et al.ll201ll ; IFridriksson et all 120111 )? 
We therefore consider this the most likely interpretation of 
the observations. Thus, we propose that the crust of the 
neutron star in J 1748 was heated during its recent 10- week 
accretion phase a nd is currently c ooling, supporting previ- 
ous speculations (JDegenaar fc Wiinandsll2011al lbh. 



1 We note that if residual accretion occurs in quiescence, varia- 
tions in the thermal emission could be due by a changing emission 
area (a hotspot at the magnetic pole) and may not necessarily im- 
ply a changing surface temperature. Given the data statistics, it 
would not be possible to make this distinction. 



4 N. Degenaar E. F. Brown & R. Wijnands 



3.1 Crust cooling simulations of IGR J17480-2446 

We employed t he thermal evolution code of 

iBrown fc Cummind (|2009T l. using all the basic physics 
ingredients and modeling approach described in that work, 
to simulate neutron star crust cooling curves for J1748 
(Fig. [3]). We assumed that the outburst had a duration 
of t b = 0.17 yr and that mass was accreted at a rate of 
M = 2.0 x 10 17 g s -1 , as infe rred from X-ray observations 
l|Degenaar fc Wijnands! l2011af ). Furthermore, the tempera- 
ture of the neutron star core was fixed at T = 7.0 x 10 7 K, 
which matches the inferred va lue of the quiescent base level 
|Pegenaar fc Wiinandsll2011bl ). 

Using standard physics input, the crust temperature is 
too low and the slope of the cooling curve too flat to match 
the observations (dotted curve in Fig. [3} • Motivated by re- 
cent theoretical conjectures (see Sec. [TJ, we also performed 
simulations with a heat source placed at shallow depth inside 
the neutron star and varied its strength, Qcxtra- We position 
this extra heat at a column depth of Pj g — 4.5X10 11 gem -2 , 
which roughly corresponds to a density p~4x 10 8 g cm" 3 
or ab out 10 m inside the neutron star IjBrown fc Cummina 
I2009T ). This is roughly the depth at which superbursts are 
thought to ignite. As can be seen in Fig. [3] adding extra 
heat raises the temperature and steepens the slope of the 
crust cooling curve (dashed and solid lines). 

For the assumed M ~ 2 x 10 17 g s" 1 
IjDegenaar fc Wijnands! l2011bi ), we find that a heat 
source of Qextra — 1-0 MeV per accreted nucleon can match 
the observational data of J1748 (solid curve in Fig. [3}. 
This is substantial compared to the total heat release of 
~ 2.0 MeV per accreted nucleon th at is accounted for 
by current nuclear he ating models IjGupta et al.l 120071 ; 
lHaensel fc ZdunikH2008l ). However, the input value of the 
mass-accretion rate has a notable effect on the resulting 
crust cooling curve and hence on the strength of the inferred 
additional heat source. If we increase M by a factor of two 
in our model calculations, which is not implausible given 
the uncertain ties in determining th is parameter from X-ray 
observations (jin't Zand et al.l 120071 ). the magnitude of the 
required extra heat decreases to Qoxtra — 0.5 MeV per 
accreted nucleon. Nevertheless, even if J1748 accreted at 
the Eddington limit, additional heat is required to explain 
the observations as cooling of the accretion-heated crust. 

As discussed bv lBrown fc Cummina (|2009|), heat is pri- 
marily conducted in the neutron star crust by electrons, 
and hence the thermal conductivity is determined by the 
rate of scatterings between electrons and ions, and between 
electrons and lattice impurities. The level of impurities in 
the crust is parametrized by the distribution of ion charges. 
For the calculations presented in Fig. [3j this was set to 
(Z 2 ) — (Z) 2 — 5, which is consistent with the values in- 
ferred from m odeling the crust cooling curves of two other 
neutron stars (IBrown fc Cumrning 2009^ and molecular dy- 
namics simulations l| Horowitz et al . 2009). We explored sim- 
ulations with different levels of impurities, but found that 
this does not affect the resulting crust cooling curve of J1748. 
This is because, unlike systems with accretion episodes that 
last for years and heat the entire crust, for short outbursts 
the heating primarily occurs at low densities, where the heat 
conduction is controlled by electron-ion scattering. 

We calculated crust cooling curves varying the depth of 



Table 1. Results from fitting the Chandra spectral data. 



Date 



kT° 



2003-07-13/14 
2009-07-14/15 

2011-02-17 
2011-04-29/30 



72.0 ±3.3 

100.9 ±3.3 

92.1 ±3.1 



1.64 ±0.30 0.59 ±0.11 



6.33 ±0.83 
4.39 ±0.60 



2.29 ±0.30 
1.59 ±0.21 



Note. kT°° is given in eV, -F ml abs in 10 13 erg cm 2 s x 
(0.5-10 keV) and L bol in 10 33 erg s" 1 (0.01-100 keV). 

the extra heat source. We found that the exact depth pri- 
marily affects the cooling rate at early times, within a few 
tens of days after the end of the accretion outburst. J1748 
could not be observed with Chandra on such a time scale, 
because the source was too close to the Sun at that time. 
Thus, the location of the extra heat source does not affect 
the crust cooling curve at the time probed by the observa- 
tions of J1748. This also implies that our Chandra obser- 
vations cannot firmly constrain the depth of the additional 
heat release, although we can obtain some limits. 

Assuming that the crust cooling curve of J 1748 indeed 
follows a power law decay with an index of —0.1 (see above), 
we estimate that there is an inward-dire cted heat flux in the 
crust of F ~ 3 x 10 21 erg cm -2 s~ x (jBrown fc Cummina 
I2009I . eq. 12). At any given depth inside the neutron star, 
it takes a certain characteristic time for heat to diffuse 
to the surface. Therefore, at a given time there will be a 
certain depth at which the layer above has already ther- 
mally relaxed, while d eeper layers are still hot and have 
not yet started to cool (|Brown fc Cumrn ing 2009). The first 
Chandra observation of J1748 was obtained ~ 50 days af- 
ter the end of the accretion outburst: this is the charac- 
teristic thermal time scale for layers located at a column 
depth of roughly P/g ~ 10 14 g cm -2 , which corresponds 
to a matt er density of p ~ 3 x 10 1 g cm -3 or a depth of 
~ 150 m (IBrown fc Cummina I2009T ). Since the heat flux is 
directed inwards, there must thus be a source of heat exte- 
rior to the depth corresponding to the thermal time scale 
of 50 days. This implies that there is a heat deposit at less 
than ~ 150 m inside the neutron star (Brown fc Cummina 
120091 ) . This is much shallower than the depth a t which most 
crust nuclear reactions are thought to occur (IGupta et al.l 
120071 ; lHaensel fc Zdunikll200Sl ; IHorowitz et al.ll2008J) . 



3.2 Repercussions 

The X-ray observations of J 1748 obtained after its recent 
10-week accretion outburst can be explained as cooling of 
the accretion-heated neutron star crust. However, regard- 
less of the exact value of M, this requires a substantial 
source of heat located at shallow depth inside the neu- 
tron star that is not accounted for by existing nuclear 



tron star tnat is n ot a ccounted tor py existing nuclear 
phys i cs calculations J Gupta et al.l 120071 ; lHaensel fc Zdunikl 
120081 ; IHorowitz et al.ll2008l ). A shallow heat source has also 
been proven necessary to explain the crust cooling curve of 
at least one of the quasi-per sistent LMXBs (MXB 1659- 
29; IBrown fc Cummina 120 09). The presence of substantial 
heat sources in the outer crustal layers is important for 
the ignition of superbursts: the extra heat release raises 
the crust temperature and can thus compensate for the 
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Figure 3. Evolution of the neutron star temperature of IGR 
J17480-2446, as inferred from analysis of Chandra spectral data. 
The grey shaded area marks the quiescent base level observed 
prior to the accretion activity of 2010 October-December. The 
three curves indicate neutron star thermal evolution simulations 
for different amounts of heat production in shallow crustal layers. 



effect of a high thermal conductivity l|Gupta et al.l 120071 : 
iBrown fe Cummindl2009h . 

Our work on J 1748 suggests that the crust of a neu- 
tron star can become substantially heated during an accre- 
tion episode of only a few weeks, and that it is technically 
feasible to detect the subsequent crust cooling with exist- 
ing, sensitive X-ray instruments. LMXBs in which the neu- 
tron star accretes for weeks or months are much more com- 
mon than those that undergo extended accretion episodes of 
years. Future crust cooling studies can therefore draw from 
a larger sample of transiently accreting neutron stars. The 
most promising targets are neutron stars that have rela- 
tively low core temperatures, but exhibit bright accretion 
episodes so that a substantial amount of matter is con- 
sumed and a considerable tempe rature gradient can develop 
in the crust (|Brown et al.lll998T) . Transient Be/X-ray bina- 
ries, which harbor strongly magnetized neutron stars feeding 
of the matter expelled by their massive Be co mpanions, can 
potentially also be included in future studies (|Brown et al.l 
ll99Sft . 

Studying a larger sample of crust cooling curves allows 
to investigate whether the presence of strong, shallow heat 
sources is a common feature amongst transiently accreting 
neutron stars. In particular, observations performed within 
~ 1 month after the end of an accretion outburst can fur- 
ther constrain the magnitu de and depth of such heat sources 
IjBrown fe Cumming 120091 ). The next challenge is then to 
identify the nature of this extra heat release. 

Our crust cooling simulations suggest that IGR J 17480- 
2446 continues to cool within the next year. New Chandra 
observations are planned to confirm this and to put further 
constraints on the shape of the crust cooling curve. 
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